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We show that the complicated band reconstruction near Fermi surfaces in the magnetically ordered 
state of iron-pnictides observed by angle-resolved photoemission spectroscopies (ARPES) can be 
understood in a meanfield level if the intrinsic unit cell doubling due to As atoms is properly 
considered as shown in the recently constructed S4 microscopic effective model. The (0,7r) or (7r,0) 
col-linear antiferromagnetic (C-AFM) order does not open gaps between two points at Fermi surfaces 
linked by the ordered wave vector but forces a band reconstruction involving four points in unfolded 
Brillouin zone (BZ) and gives rise to small pockets or hot spots. The S4 symmetry naturally chooses 
a staggered orbital order over a ferro-orbital order to coexist with the C-AFM order. These results 
strongly suggest that the kinematics based on the S4 symmetry captures the essential low energy 
physics of iron-based superconductors. 

PACS numbers: xxxx 



Introduction: The parent compounds of all Iron- 
pnictide high temperature superconductor p]-[3| are char- 
acterized by a unique C-AFM order 0-0 • In the past 
several years, although the majority of theoretical stud- 
ies on these compounds have successfully obtained such 
a magnetic state from both itinerant and local moment 
models [Ml, the band reconstruction in the C-AFM 
state observed by ARPES jl5Ml8| remains a puzzle. Un- 
like in a conventional spin density wave (SDW) state 
where gaps at any two points at Fermi surfaces which 
are connected by the ordered SDW wave vector should 
be opened generically, the band structure in the C-AFM 
state appears to be reconstructed by the development of 
the magnetic order. While optical measurements suggest 
a significant portion of Fermi surfaces are gapped[19l-l22|, 
ARPES observes that small pockets or hot spots develop 
around Fermi surfaces. It is fair to say that the sce- 
nario of the gap opening due to a magnetic order in a 
multi-orbital system, such as iron-pnictides, is complex. 
For example, it has been argued that the entire Fermi 
surfaces can be gaped out [23]. Nevertheless, the surpris- 
ing band reconstruction is never clearly understood until 



In this paper, we demonstrate that in the recently con- 
structed effective two-orbital model with the S4 symme- 
try, which is the symmetry of the trilayer FeAs struc- 
ture, the band reconstruction becomes a natural conse- 
quence from the intrinsic unit cell doubling due to As 
atoms [24]. The spirit of the S4 model is that the kine- 
matics can be divided into weakly coupled two orbital 
models controlled separately by two S4 iso-spin compo- 
nents as shown in FigUfa). For each S4 iso-spin compo- 
nent, a unit cell includes two irons, namely, the unit cell 
is doubled, and the doubling results in the formation of 
a pair of bands which are responsible for a hole pocket at 
T and an electron pocket at M. Although this kinemat- 
ics is hidden in various multi-orbital models constructed 
for iron-based superconductors as shown in [24]. Here we 
show that it is the extraction of such a minimum effec- 
tive model allows us to naturally capture the band recon- 
struction features mentioned above in a simple mean- field 
Hamiltonian. The S4 symmetry also naturally imposes a 
staggered orbital order over a ferro-orbital order to co- 
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FIG. 1: (color online) (a) The lattice structure of Fe-As tri- 
layers. The different filled green/big balls indicate the top 
and bottom As layer. The bule/deep and white/light orbital 
patterns indicate two S4 iso-spin components. The effective 
hopping parameters up to 3rd nearest neighbor are marked, 
(b) The (tt, 0) C-AFM order are shown. Four Fe sublattices 
in the unit cell of C-AFM state are marked, (c) and (d) The 
Fermi surfaces and energy spectrums along high symmetric 
lines for normal paramagnetic state are shown. The hopping 
parameters are set t\ — 0.47, t[ = 0.53, £2 = 0.9, t' 2 = —0.18, 
£3 = 0.0, £3 = 0.1, t c = 0.1 and \i — —0.5 for electron under- 
doped 4.5%. 



exist with the C-AFM order. Our results suggest that 
a proper starting point in kinematics is crucial to obtain 
a correct meanfield result and the kinematics based on 
S4 symmetry due to As atoms provides the essential low 
energy physics of iron-based superconductors. 

Oriented gap opening mechanism with unit cell dou- 
bling: The effect of the unit cell doubling on the C- 
AFM state can be elucidated by considering the kine- 
matics of a single S4 iso-spin component. Deeply in a 
C-AFM state, regardless of the origin of the magnetic 
ordering, we should expect that a SDW order can cap- 
ture the basic features of itinerant bands. As shown in 
FigOJa), the kinematics of the single blue (deep) orbital 
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FIG. 2: (color online) (a) The energy bands of Hamiltonian 
Eq. (pQ) with spin-up component in (tt, 0) C-AFM state along 
r — M x and T — M y directions are shown. The bule/dashed 
lines indicate the normal state while the red/solid lines in- 
dicate the C-AFM state. The shadow regime indicates the 
parts below Fermi surfaces, (b) The shadow regime in (a) is 
zoomed in. The arrows indicate the bands in C-AFM state 
shift up or down compared with that in normal state. We set 
ti t[ 0.5, t 2 0.9, t' 2 -0.18, t 3 = t' 3 = 0.0, fi -0.69 
and Asdw — 0.15. 



has an intrinsic unit cell doubling due to the large dif- 
ference between the two next nearest neighbor hopping 
parameters, t<i and t' 2 (note: the details of the S4 model 
can be found in our previous paper |24j). 

In the Qi = (tt, 0) C-AFM state shown in Fig. 1 (b), 
at the 1 and 4 sites, the spin- up component is a ma- 
jority while at the sites 2 and 3, the spin-down compo- 
nent is a majority. A simple SDW order is thus given 
by Asdw =< Hka^ c \,a c k+Qi^ > where r] t = 1 and 
rj^ = — 1 for different spin components. Due to the in- 
trinsic unit cell doubling, for any given point k in un- 
folded BZ zone, there are four moment urns linked in the 
presence of the magnetic order, fc, k + Qi, k + Q, and 
k + Q2 where Q = (7r,7r) and Q2 = (0,7r). The mean- 
field Hamiltonian that described the (71", 0) C-AFM state 
can thus be expressed as H = ^2 k a ^ a l-L(k, cr)^^ with 
^k,a = [cfe,a,c fe+ Q 5a ,c/ c+ Q l5a ,c/ c+ Q 2 , a ] T . U(k,o) is given 



by 
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where /(fc) = 2ti s (cos fc^+cos fc 2/ )+2tid(cos fc^ — cosk y ) 
-\-At2s cos k x cos k y + 2t3 S (cos2/c :c + <zos2k y ) 
-\-2t^d{cos2k x — cos2k y ) — /i, g(k)=At2d sin fc^ sin fc^ 
and = (** + *£)/2, t id = (t 2 - *' 2 )/2 with z = 1,2,3. 
It is explicit to find that the spin-up and spin-down 
electrons play a symmetrical role in building up the 
(7r, 0) C-AFM states. Hence, we can only investigate 

«(*,t). 

Now we are in a position to show that the gap opening 
described in EqJT]is rather different from a conventional 
one in a conventional SDW state. Here the energy disper- 
sions of H(/c,t) can n °t be given by a simple analytical 
form. We provide the numerical results for spin-up case 
shown in Fig. 2. It is easy to see that the band recon- 
struction in the C-AFM state are anisotropic along the 
T — M x and Y — M y directions. Below Fermi surfaces, an 
energy splitting takes place between two hole bands at 



(f , 0) point along the Y - M x . Along the Y - M y direc- 
tion, gaps at Fermi surfaces are opened between a hole 
band and an electron band. Due to the fact that the gap 
openings are different along different directions, we call 
this behavior as oriented gap opening. As we will show 
later, this behavior leads to a development of hot spots 
or a small pockets on Fermi surfaces. 

The self- consistent mean- field results: In the above 
analysis, the magnetic order is introduced artificially. To 
show the magnetic order is generated by standard inter- 
actions in the S4 model (note, a recent quantum monte 
carlo calculation has already shown the C-AFM correla- 
tion is the leading magnetic correlation of the model near 
half filling [11]) as well as to produce more realistic results 
that can be closely compared to ARPES observations, we 
apply a self-consistent meanfield calculation to consider 
the full S4 model with the full interaction terms. 

Including all possible electron-electron interaction 
terms, the general Hamiltonian of the S4 model is given 
by Hi, 

H = Ho + H int (2) 
where Hq is the kinematics part given in[24| and 



H^t = U y^(nL+nLi + nL^n d 
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Here, we use c and d to identify the two S4 iso-spin com- 
ponents. h c aia = cl ia c aia and h d ai(T = dl ia d ai(T . a la- 
bels the two sublattices of the iron square lattice with 
a = 1, 2. U and U' are the direct intra- and inter-orbital 
Coulomb repulsions, respectively. Jh is the Hund cou- 
pling satisfying U' = U — 2J#. The last term in Hi n t 
describes the spin-flip term of intro-orbital exchange and 
the pair hopping. Here we neglect the pair hopping term. 

To obtain the possible ground-state in the undoped 
and under doped cases, we apply a mean- field decoupling 
to the model in Eq.(|3j). H in t can be simply decoupled by 
defining the mean- field values of the diagonal operators, 
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If we assume that possible ordered states are only allowed 
to break translation symmetry up to four irons and let 
a = 1,2,3,4 label sites in the four different sublattices 
respectively as shown in FigQJb), the mean- field Hamil- 
tonian of Eq.([2]) in the reciprocal space can be written 
as 
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where C represents non-operator terms. 

The above mean-field Hamiltonian can be numerically 
solved with a standard self-consistent procedure by mini- 
mizing free energy. During the iterative procedure of the 
self-consistent calculation, we enforced the average elec- 
tron number per site at each step remains a constant, 
namely, keeps the density constant. We define the mag- 
netization m a and polarization A 77/q/ as 

m a = « T - n c al ) + (n d at - n d ai ) (7) 
An a = « f + n c al ) - « T + n d al ) (8) 

All different kinds of standard magnetic orders (such 
as C-AFM order, ferromagnetic order and antiferromag- 
netic order) and orbital orders (such as ferro-orbital or- 
der, staggered-orbital order and stripe orbital order) are 
included by above meanfield ansatz. 

The zero-temperature phase diagram for the undoped 
case was obtained and shown in Fig. 3 (a). The C- 
AFM order was found to be stable in a broad and re- 
alistic range of interaction strength. Furthermore, the 
C-AFM orders exhibit differently depending on U and 
Jh/U. When Jh/U is large in a moderate value of [/, no 
orbital-order was found. When Jh/U is small enough, a 
small staggered-orbital order coexists with the magnetic 
order. When U is large enough, as expected, the metal- 
lic C-AFM order becomes insulating. It is important to 
note that for a coexisting state with both C-AFM order 
and staggered orbital order, the corresponding magneti- 
zations m a at the sites labeled in FigJTJb) are m i/4 — m 
and m 2 /3 = — m and the orbital orders are given by 
An; = m (i = 1,2,3,4). In Fig. 3 (b) (c) and (d), 
we plotted m and m as functions of U, Jh/U and tem- 
perature T. From Fig. 3 (b), for the two undoped and 
4.5% electron doped cases studied here, m jumps into a 
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FIG. 4: (color online) The Fermi surfaces in C-AFM states 
are shown for some different interaction strength and doping, 
(a) the undoped case with [7 = 4 and Jh / U = 0.15. (b) 
the undoped case with U = 3.5 and Jh/U = 0.15. (c) the 
4.5% electron-doped case with [7 = 4 and Jh/U = 0.25. (d) 
(e) and (f) are the energy bands along F — M x and F — M y 
directions, respectively. 




FIG. 3: (color online) (a) The zero-temperature mean- field 
phase diagram of the S4 effective model in the Jh/U vs U 
plane in undoped case. We denote 'PM' for the paramagnetic 
phase, 'C-AFM V for the metallic C-AFM phase, 'C-AFM 2' 
for the metallic C-AFM phase with a staggered orbital order 
and 'C-AFM 3' for the insulating C-AFM phase, (b), (c) and 
(d) show the C-AFM magnetization m and stagger orbital 
polarization m as the function of U, Jh/U and T. The red 
color denotes the undoped case while the blue color denotes 
the 4.5% electron-doped case. The parameters are Jh/U = 
0.15 for undoped case and 0.25 for 4.5% electron doped case 
in (b), U — 4 for both cases in (c), and Jh/U with the same 
values as (b) and U — 4 in (d). The hoping parameters are 
indicated in the caption of Fig. 1. 



finite value at U ~ 3.2 and ~ 3 respectively, which means 
the transition from the PM to the C-AFM is a first order 
phase transition. From Fig. 3(c), the large Jh/U en- 
hances the C-AFM order m and suppresses the staggered 
orbital order m in both undoped and underdoped cases. 
From Fig. 3 (d), increasing temperature suppresses both 
C-AFM order m and the staggered orbital order m G , and 
so does the doping. We also find that increasing doping 
can continually reduce the C-AFM order m at low tem- 
perature. All the results are qualitatively consistent with 
the experimental observations [25j. 

The result that a staggered orbital order is favored over 
a ferro-orbital order stems from the S4 symmetry. As 
shown in FigJTJa), eacn S4 iso-spin component is formed 
by two different types of orbitals in the two different sub- 
lattices of the iron square lattice. A staggered orbital 
order corresponds to the S4 symmetry breaking. Many 
theoretical proposals based on models with one-iron per 
unit cell emphasize the ferro-orbital order [26l-[28j. There- 
fore, in principle, the result can be used as a test for the 
S4 symmetry. Unfortunately, the staggered orbital order 
is rather small in our meanfield results. 

The formation of hot spots or small pockets in the C- 
AFM state: To obtain further insights into the charac- 
ters of the C-AFM ordered states, we investigated the 
Fermi surface construction which can be directly ob- 
served by APRES experiments. The reconstructed Fermi 
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FIG. 5: (color online) The band structure around (f ,0) cor- 
responding to Fig 4. (a) (d). The Dirac cone structures are 
marked with the red-dashed lines. 

surfaces and bands of our model with different inter- 
action strengths and doping cases are shown in Fig. 
4. The interaction strength is chosen in the moder- 
ate strength range since iron-pnicitides are most likely 
in intermediate coupling region. The Fermi surfaces in 
Fig. 4 (a), (b) and (c) show a close match to ARPES 
results measured in BaFe2As2[16], NaFeAs[l7l fl8j and 
Ba(Fei_^Co^)9As9fl5j]. respectively. Furthermore, the 
magnetization value in Fig. 4 (a) with m = 0.48 for 
BaFe2As2 is larger than the one in Fig. 4 (b) with 
m = 0.33 for NaFeAs, which is also quantitatively con- 
sistent with experiments (25|. 

The remarkable character of reconstructed Fermi sur- 
faces in BaFe2As2 is the Dirac cone structures around 
the (j ,0) and (^,0). Especially, the energy spectrum 
around the Dirac cone structures in Fig. 4 (d) are sim- 
ilar to the experimental fitting results [16]. From Fig. 4 
(d), we can find that the Dirac cone structure around 
(f,0) point is formed by a hole-like band from one S4 
iso-spin component and a folded electron-like band from 
the other S4 iso-spin component while the Dirac cone 
structure around (^,0) point is formed by exchanging 
the contributions between two S4 iso-spin components. 



In Fig. 5, we plotted the spectra around the (f , 0) point, 
where the Dirac cone structures are resolved with the 
red-dashed line. The Dirac cone structures were pro- 
posed as the universal character and robust with inter- 
action strength and doping [23, |29|. Our results clearly 
show that the feature highly depends on quantitative pa- 
rameter settings in a model and in general should be 
material-dependent . 

Discussion and conclusion: Several LDA calculations 
have also noticed that some properties have to be un- 
derstood in the base of the 2 Fe unit cell[30|, |3l|. From 
above results, we clearly see this importance reflected in 
the starting point of kinematics in analytic approaches. 
In most previous studies, in particular, analytic studies, 
calculations are performed in models with one iron per 
unit cell. In principle, the S4 symmetry kinematics is hid- 
den in these models and can be revealed by performing a 
unitary transformation. However, when interactions are 
concerned and exact solutions are impossible, the search 
of meanfield order parameters highly depend on the base 
set in kinematics, which is the key reason why an explicit 
S4 construction becomes crucial. The result also suggests 
that many physical quantities in a symmetry broken state 
require new investigation. 

It is also important to note that our calculation is per- 
formed in a tetragonal lattice. It is known that a lattice 
distortion or a tetragonal to orthorhombic lattice transi- 
tion exists at a temperature higher than or equal to the 
magnetic transition temperature [25]. In general, the S4 
symmetry is broken in the orthorhombic lattice. Our cal- 
culation does not include this aspect since the model is 
purely electronic and does not include an electron-lattice 
coupling. The lattice distortion is likely driven by pure 
electronic nematismJH, 1321434) so that the qualitative re- 
sults here will remain unchanged. Nevertheless, a full 
consideration of this aspect is needed in future. 

In summary, we show that the S4 microscopic effective 
model characterized by a doubling of iron unit cell ex- 
plains the puzzled band reconstruction in the magnetic 
state of iron-pnictides. The results provide a strong sup- 
port to the kinematics described by the S4 symmetry. 
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